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ABSTRACT 

S'-matrices associated to the vector representations of the quantum groups for the 
classical Lie algebras are constructed. For the a m _i and c m algebras the complete S- 
matrix is found by an application of the bootstrap equations. It is shown that the simplest 
form for the ^-matrix which generalizes that of the Gross-Neveu model is not consistent 
for the non-simply-laced algebras due to the existence of unexplained singularities on the 
physical strip. However, a form which generalizes the S'-matrix of the principal chiral 
model is shown to be consistent via an argument which uses a novel application of the 
Coleman-Thun mechanism. The analysis also gives a correct description of the analytic 
structure of the S'-matrix of the principle chiral model for c m . 
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1. Introduction 



The complete non-perturbative solution of interacting quantum field theories seems 
like a hopeless fantasy. In the vast space of two-dimensional quantum field theories, how- 
ever, there is a small set of theories attracting a large amount of interest due to the fact 
that they are integrable and hence to a certain degree solvable. At the moment the under- 
standing of such theories is mostly at the level of on-shell physics encoded in the scattering 
matrix. Even though the S-matrix seems to be rather simple, since it factorizes and is 
specified completely by the two-body S-matrix, finding the complete S-matrix satisfying 
all the required properties is a surprisingly difficult task. Also it is worth pointing out that 
a given S-matrix may always be multiplied by CDD factors to give an equivalent S-matrix. 

Let us summarize the situation so far. First of all, there exist a set of minimal purely 
elastic S-matrices (and hence trivial solutions of the Yang-Baxter equation) which are 
related to the simply-laced affine Lie algebras. There are particular CDD factors, which 
depend on a coupling constant which make these S-matrices equal to the conjectured 
S-matrices of the corresponding affine Toda field theory (the connection with the field 
theory being established in perturbation theory) [1]. There are also a set of S-matrices 
which describe the non-simply-laced affine Toda field theories, however, the situation is 
rather different from the simply-laced theories in that the S'-matrix does not factor into a 
'minimal' piece independent of the coupling and furthermore the mass ratios depend on 
the coupling constant [2]. 

We now turn to the S-matrices that are non-trivial solutions of the Yang-Baxter 
equation, so that the particle states carry some internal quantum numbers. One has both 
the rational and trigonometric solutions of the Yang-Baxter equation to hand, the latter 
depending on a coupling constant. S-matrices corresponding to rational solutions have 
been considered in [3], and describe the principal chiral model for a classical Lie algebra 
and the Gross- Neveu model for the simply-laced classical Lie algebras. It is worth pointing 
out that a Gross-Neveu type S-matrix for the non-simply-laced algebras can be written 
down, but the Ansatz fails due to the existence of unexplained singularities on the physical 
strip. This point highlights the delicate nature of these S-matrices; it is simply not good 
enough to write down an Ansatz for the S-matrix elements on some elementary particles, 
hoping that the bootstrap will close on a conjectured set of particles; in order to claim a 
consistent S-matrix it is necessary to account for all singularities on the physical strip in 
terms of bound states or anomalous thresholds via the Coleman- Thun mechanism [4] . 

Trigonometric solutions of the Yang-Baxter equation are more general since they de- 
pend on a coupling constant and the rational solutions can be obtained from them in a 
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certain limit. S matrices constructed from trigonometric solutions have been the subject of 
much speculation but very little detailed analysis, except for the series associated to a m _i 
[5,6] and in particular for a\ which gives the soliton S'-matrix of the sine-Gordon theory 
[7] . In these cases one can prove that the bootstrap closes on a given set of particles, when 
the coupling constant satisfies a certain inequality [5]. The problem encountered for other 
algebras stems from the difficulties in solving the bootstrap equations. However, we shall 
show how this can be achieved for c m (in addition to a m _i which has been considered 
previously). 

One of the principle reasons for constructing the trigonometric S'-matrices is that they 
are thought to describe the integrable perturbations of certain conformal field theories, 
typically displaying VF-algebra type symmetries [6,8]. 

2. Trigonometric solutions of the Yang-Baxter equation 

In this section we describe the trigonometric solutions of the Yang-Baxter Equation 
(YBE). They are associated to certain deformations of the universal enveloping algebra 
of a Lie algebra known as a quantum group [9,10]. The solutions can be thought of as 
intertwiners between tensor products of representations of the algebra: 

R(u) : V^V U ^V V ®V^ (2.1) 

where u is the (additive) spectral parameter. Such an '.R-matrix' has a spectral decompo- 
sition 

£(u) = J]p A (u)P A , (2.2) 

A 

where Pa is a quantum group invariant homomorphism <g) V v — > V v <S> with the 
property that Pa|v a , 7^ if and only if A' = A. 1 If \i = v then Pa is a projection. 

In the following we shall use both the language of spectral decompositions and the 
interaction-round-a-face (IRF) picture when writing down solutions of the YBE. When 
using the latter we shall denote R(u) as W(u). To start with we consider the solutions 
associated to the vector representation of the algebra. In what follows we shall consider 
all the classical Lie algebras: a m _i, 6 m , c m and d m . The set of weights E of the vector 
representations are 2 

£ = {ei-(eiH h e m )/m, . . . , e m - (ei H he m )/m}, for a m _i, 

E = {0, ±ei,...,±e m }, for b m , E = {±e 1 , . . . , ±e m } , forc m ,d m , 

1 We are assuming that every irreducible component V\ C V M <£> V„ has multiplicity one. 

2 In our notation the long roots of b m have length 2 while for c m the single long root has length 
4. 



2 



where the e^'s are a set of ortho normal vectors. 

The solution of the YBE is labelled by four weights of the algebra: 



W 



a 



d 



a | , a,b,c,d G A*, 



(2.4) 



with the property that W is only non-zero if c — a, d — c, b — a and d — b are G E. 

For completeness we now write down the solutions following [11] (see also the review 
[12]). In the following uj is a constant which is related to the deformation parameter of the 
quantum group. For convenience we introduce for a G A* 



= uj(a + p) • p, for |i G S / 0, 
cxo = — uj/2, a^ v = — a v , a^- v = + a v , 



(2.5) 



where p is the sum of the fundamental weights of the algebra. 3 Also we define 



[x] = sinx, A = tgu/2, 



(2.6) 



where g is the dual Coxeter number of the algebra and t is the (length) 2 /2 of the longest 
root. 4 

For a m _i the solution is 



W 
W 

W 



u 



a a + p 

a + p a + 2/i 

a a + p 

a + p a + p + v 

a a + v 

a + p a + p + v 



u 



[u - \u]/[u], 

= [a^u + Aw]/[a^] 



(2.7) 



u = 



1 /2 



3 These are the vectors Ui with Ui ■ oij = {a 2 j/2)5ij where the otj are the simple roots. 

4 The values (g, t) are (m, 1), (2m — 1,1), (m + 1, 2) and (2m — 2, 1) for a m -i, b m , c m and d m , 
respectively. 
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where fj,, v e £ and v. For the algebras 6 m , c m and <i m the solutions are 



W 
W 
W 



a a + u 
a + u a + 2u 



t n = [X ~ X *r Xu] , for^O, 



a a + u 

a + u a + u + v 

a a + v 

a + u a + u + v 



u 



u = 



[X][u] 
[A - Xu] [a^v + Xu] 

[a][<v] 

[A - Xu] [Xu] 



, for n 7^ ±i/, 



[A][o;] 



x 



w 



a a + v 
a + u a 



u 



[a^ + ui\ [a^ - ui\ 

[Xu] [a M _„ + uj - X + Xu] 
[A] [a M _„ + uj] 



1/2 

, for u 7^ ±z/, 



(2.8) 



[A - Att]^-^ + uj + Xu] , 
+ V FTTTI^ —i ' for/Lt^O, 



W 



a a 
a a 



u 



[X][a^- U + uj 
[X + Xu][2X - Xu] [Xu] [X - Xu] 



where u, v G £ and 



[A] [2 A] [A] [2 A] 

[a MK + uj 



_ sja^ + uj) 
~° 8 ( a ) 11 \ a 

T _ K + w / 2 - 2A ] ^ 
Ja — — n — i rrToi — ( - Ta -- 



, for u ^ 0, G ao = 1, 



(2.9) 



K + w/2] 



In the above a is —1 for c m otherwise being 1. The function s(x) = [tx] for b m and c m , 
and s(x) = 1 for d m . 

The solutions satisfy a set of conditions, in addition to the YBE, which will be im- 
portant for the construction of an S- matrix (our terminology follows that of [12]). 
(i) The standard initial condition. 



W 



a b 
c d 



(ii) The unitarity condition. 



a e 
c d 



u W 







a b 

e d 



(2.10) 



u = g(u)5 bci 



(2.11) 



where 



. , [uj - Xu] [uj + Xu] , 

Q{u) = rT5 , for a m -i, 

[uj] 2 

[X — Xu] [uj — Xu] [X + Xu] [uj + Aw] 

= WW 



(2.12) 



, for b m , Cm and d r 



(iii) Crossing symmetry (b m , c m and d m only). 



W 

where 



a b 
c d 



u = W 



c a 
d b 



G a = e(a) JJ s(a;) j] [en - aj][oi + aj], (2.14) 

k=l l<i<j<m 

where a = Y^L\ a i e i defines the a, and e(o) is a sign factor chosen so that e(a+fi)/e(a) = a. 
G a is related to G a ^ in (2.9) by G a ^ = G a+fl /G a . 

For a m _i the representation is complex and so the solution of the YBE does not 
satisfy a crossing symmetry relation involving only the vector representation. 

The above solutions correspond to the unrestricted solutions of the YBE where the 
variables {a, 6, c, d} are any weights of the algebra. We shall be primarily interested in the 
restricted models which are obtained for the particular values 

»=WTFy (2 ' 15) 

and the weights are restricted to lie in the set of integrable weights of the affine algebra at 
level k projected onto the weight lattice of the finite algebra. In practice, this means that 
we restrict the allowed weights to the dominant weights satisfying 

a-6<k, (2.16) 

where 9 is the highest root of the algebra. We denote the set of weights satisfying this 
condition as A*(k). 



3. The vector-vector scattering matrices 



The solutions of the YBE equation that we wrote down in the last section naturally 
lead to S'-matrices for a set of kinks. We denote a kink state by K ab {9), where a and b are 
two vacua of the theory and 6 is the rapidity of the kink. In an integrable field theory we 
need only consider the S'-matrix for the process 

K ac (0i) + K cd (02) - K ah {9 2 ) + K bd (0i) (3.1) 

since all the other S'-matrix elements are determined in terms of these. The idea is to 
find an S'-matrix for a theory associated to a classical Lie algebra with vacua in one-to- 
one correspondence with A*(k) (or A* in the unrestricted model) and hence kinks with a 
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topological charge being a weight of the vector representation of the algebra. The S-matrix 
of the process (3.1) is related to the solution of the YBE as 

5(« ^^(uW^ ^|/ M )(||) /< " ,/2 . (3.2) 

for some scalar function Y(u), where u = (9\ — 62)/^ is the rapidity difference of the 
incoming kinks. The other scalar factor is included in order to simplify the implementation 
of crossing symmetry. The scalar function Y(u) must be chosen so that the S'-matrix 
axioms are satisfied, subject to minimality, meaning that we search for Y(u) which ensures 
the axioms are fulfilled with the minimum number of poles and zeros on the physical strip 
(the region < Re(-u) < 1). 

The unitarity constraint which follows from the hypothesis that the space of states 
of the theory is complete can be satisfied by virtue of the completeness relation (2.11) if 
f(u) = cu, for some constant c, and 

Y(u)Y(-u) = 1/q (/(«)). (3.3) 

Crossing symmetry of the S'-matrix relates the process (3.1) with the the process K c d + 
Kdb — ► K ca + K a b, where K is the charge conjugate kink. Notice, that for the b m , c m and 
d m algebras if a — b G E then b — a G S, due to the fact that the representations are real, 
therefore we may take K ab = K ah . In these cases crossing symmetry requires 



S 



a b 
c d 



u)=S 



c a 
d b 



1 - u ), (3.4) 



which may be satisfied by virtue of (2.13) if f(u) = u and 

Y{l-u) = Y{u). (3.5) 

The situation for a m _i is somewhat different since the vector representation is not 
conjugate to itself. However, this case has been dealt with elsewhere [5,6]. One finds 
that a crossing symmetry relation can be satisfied when kinks in the conjugate vector 
representation are included in the spectrum and the minimal solution for Y(u) is 

dx 2 sinh(mwa;/2) 



Y m ,k(u) = exp {/ 



x 



x sinh[(/c + m)x]smh(mx) (3 6) 

(cosh(/ca;)cosh(ma;w/2) — cosh[(m + k — 2)x]cosh[mx(u/2 — 1)]) j, 



where k is a parameter related to uj as in (2.15). Y m ^{u) can also be expressed in terms of 
products of Gamma functions [6]. The function Y m ^{u) satisfies two important identities: 



sin 2 uj 



Y m , k (u)Y mjk (-u) = — — . — , (3.7) 

sin (a; + Xu) sm(u — Xu) 
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and 



y m , fc (i - u)y m ,fc(i + u) 



sin 2 u; 



(3.8) 



sin(A — Aw.) sin(A + Aw) 

where as before u = ir/(m + k) and A = rruo/2. 

Using these two relations we can now write down S'-matrices for the other algebras 
by noting that a solution of (3.3) and (3.5) is 



Y(u) = Y t g,tk(u)Y t g,tk(l ~ i*)-?— ~> for 6 m ,c m ,rf m . 
y y ' sin a; 



(3.9) 



Summarizing we have 



u 



= W 



u 



X 



GgGd 
GhG r 



Y m ,k{u), 

Y t g,tk{u)Y tg , tk {l - u) 



u/2 



sin A 



for a m _i, 
for b m , c m , d m . 



(3.10) 



It is straightforward to show using (2.10) and the fact that Y tg ^{Q) = 1 and Y tg ^{l) = 
sincu/ sin A that 

a b 



S 



c d 



= 5 bc . 



(3.11) 



Next we turn to an investigation of the analytic structure structure of the above S'- 
matrices. The factor Y tg ^{u) has no poles or zeros on the physical strip which means that 
there are no bound state resonances. So the S-matrix defines a complete theory for the 
cases bm, c m and d m with kinks associated to the vector representation of the algebra. For 
the the S-matrix does not make sense on its own because the requirement of 

crossing symmetry means that kinks in the conjugate representation should be included in 
the spectrum. 

We now consider the case of the a m -i theory in more detail. In order to satisfy 
crossing symmetry the conjugate vector representation must be included in the spectrum. 
This is achieved in a dynamical way by a generalization of the original Ansatz (3.10): 



d 



u = X(u)S 



d 



u 



(3.12) 



where the additional CDD factor X(u) does not upset the unitarity or crossing symmetry 
relations but provides an additional pole on the physical strip. In order to motivate the 
form of the scalar factor consider the spectral decomposition of the .R-matrix (see appendix 
A): 

R {u) = 5in( l- A "» P 2 „, + Sin( l +A "V 2 , (3.13) 



sincu 



sin a; 
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where f > 2 LOl and P W2 are the quantum group invariant projectors onto the representations 
with highest weights 2ui and U2 which appear in the tensor product of two vector represen- 
tations. The idea is that X{u) should have a simple pole at u = uj/\ = 2/m at the place 
where the .R-matrix projects onto the second fundamental representation: R(u) oc P W2 . 
This determines the form of the scalar factor: 



X ( u ) = ■ )L Ty ( 3 - 14 ) 

As a consequence a particle transforming in the second fundamental representation in 
included in the spectrum and the S-matrix elements of this new state can then be found 
using the bootstrap equations. The process is then repeated: the new S'-matrix has simple 
poles on the physical strip which are interpreted in terms of states propagating in the direct 
or crossed channels. For the unrestricted models it was shown in [5] that the procedure 
terminates on a finite set of particles corresponding to all the fundamental representations 
of a m -i (in particular the conjugate vector representation) if the coupling constant satisfies 
the inequality uj < 2n/m. Notice that the restricted models, for which uj = n/(m + k), 
all lie in this region. The most important comment to make about this procedure is that 
the positions of the poles are dictated by the fusion structure of the solutions of the YBE, 
which in turn determines the ratios of the masses to be 

m a oc sin ( — ) , a = 1, 2, . . . , m — 1. (3.15) 
V m / 

Here the particle with mass m a transforms in the a th fundamental representation. This 
is the mass spectrum of the S'-matrix of the a m _i affine To da field theory. This is not 
surprising because the extra CDD factor in (3.12) is the minimal 5 S- matrix of the afore- 
mentioned theory for the particle associated to the vector representation. In fact the poles 
on the physical strip of the non-diagonal theory are completely determined by the minimal 
To da S-matrix since the rest of the S-matrix has no poles or zeros on the physical strip. 
It is important to point out that the minimal Toda scalar S-matrix has double poles on 
the physical strip, but these are explained in terms of the anomalous thresholds produced 
by 'on-shell' diagrams [1]. We will discuss the full solution of the bootstrap equations in 
the next section. 

We can now try to construct S-matrices for the other classical Lie algebras along the 
same lines as for the a m _i algebra. The first thing to consider is the spectral decomposition 
of the solution of the Yang-Baxter equation for vector-vector scattering. It is known that 



By minimal we mean the part of the Toda S'-matrix which is independent of the coupling 
constant. 
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[9] 

~ . . sin(A — Xu) sin(u; — Aw) sin(A — Xu) sin(u; + Xu) 

R{u) = : : — r 2uJl H : : — 

sin u; sin A sin a; sin A 



( Q 1 ft \ 

sin(A + Aw) sin(u; + aXu) 

H : : — , ™0, 

sin w sin A 

where cr is = — 1 for c m and 1 otherwise. So if we wish to include the second fundamental 
representation in the spectrum we would require a simple pole at u = u/X = 2/tg; however, 
notice that for c m the residue at this point projects onto the irreducible representation, 
whereas for the other algebras the residue projects onto a reducible representation consist- 
ing of the second fundamental representation plus a trivial representation: we are forced 
to conclude that this second particle transforms in this reducible representation in these 
cases. A simple pole at this point can be introduced by including a CDD factor 

sin (zp + sin ff + f - 

XW = 2 " " 2 2 (3.17) 

™ " sin - § + =) 

which also has a simple pole corresponding to the process where the second particle is 
exchanged in the crossed-channel as demanded by crossing symmetry. The form of X(u) 
ensures that it satisfies unitarity and crossing independently of the non-diagonal part of 
the S-matrix. 

So the proposal for the S-matrix for vector-vector scattering has the form 

S(w) = X(u)S(u), (3.18) 

where S(w) is the S-matrix in (3.10) and X(u) is the CDD factor (3.14) or (3.17) which is 
the minimal S'-matrix for the vector-vector scattering of an associated affine Toda theory. 
The associated Toda theories are the following: 

a m — 1 Q>m — li b m > 4m-l' C m > > ^rn ■ (3.19) 

The associated algebra is in fact the one whose Cartan matrix is equal to the transpose 
of the Cartan matrix of the untwisted affinization of the original algebra. We emphasize 
that X{u) is not the Toda S'-matrix itself but rather the part which is independent of 
the coupling constant (for the non-simply laced algebras we refer to the naive S-matrix 
elements that one would write down for the particles with the classical mass ratios [1] as 
opposed to the actual S-matrix of the Toda theory [2]). 

On the basis of this there are serious objections to the S-matrices that we have pro- 
posed. Firstly, for b m and d m the vector representation cannot be considered as the 
elementary particle from which all the others follow as bound states because the spinor 
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representations could not produced in this way. One might think that one could define a 
consistent S-matrix which does not include the spinor representations and for which the 
vector representation is elementary, however this cannot be so because the spinor makes 
an appearance in the correct interpretation of the higher order poles. In other words such 
an S'-matrix would have some unexplained higher order poles and must consequently be 
rejected. 6 The second thought is that if one starts with the spinor representations in a 
similar way then the S-matrices that we have written down would follow by fusion, in 
other words we have only written down part of the S'-matrix. 

There is, however, a more serious objection which applies to the non-simply-laced 
theories. The objection arises from the fact that for the non-simply-laced algebras the 
minimal Toda S-matrix part of the conjecture does not by itself define a consistent S- 

(2) (2) 

matrix. The reason is that as written down the scalar S-matrices for the a^m-i an< ^ ^m+i 
Toda theories have multiple poles on the physical strip which cannot be explained by the 
spectrum of particles which follows from the classical Toda Lagrangian [1] : 



j = m 

2 '"" 1 ' "' J A 1 2 sin T :-^ T j = l,2,...,m-l, 



(2) 




(3.20) 



d ™+i : m 3 « sin K7^rrr\ » j = i,2,...,m. 



The best way to see this is to consider the a^m-i an< ^ ^m+i S-matrices as a subset of the 
minimal S-matrix of a larger Toda theory. The a^m-i S-matrix is the S-matrix for 
the subset of particles {2,4,..., 2m}, whilst the d^ +1 S-matrix is the d^ +2 for the subset 
of particles {1,2, ... ,m} (i.e. without the spinor and anti-spinor particles). The point 
now is that the sub-S-matrix has multiple poles on the physical strip which require the 
full spectrum of the larger theory in order to be explained. 

In fact the resolution of this problem for the non-simply-laced Toda theories is sur- 
prisingly subtle. What happens is that the ratios of the particle masses get renormalized 
in a non-trivial way such that the mass ratios depend on the coupling constant with the 
consequence that the analytic structure becomes modified. The new S-matrix no longer 
suffers from the pathologies of the naive one. The non-diagonal S-matrices under discus- 
sion would, on the face of it, suffer the same fate as their naive non-simply-laced Toda 
counterparts, since it is the naive Toda factor that provides the poles on the physical strip; 
however, in this case we cannot hypothesize that the particle masses get renormalized in 
some non-trivial way because the positions of the poles must match the fusing of the S fac- 
tor and this fusing structure is completely rigid: the positions of the poles are determined 
completely. 



6 It is just possible for the matrix factor of the 5-matrix to provide some judicious zeros, but 
this is not the case. 
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If the S-matrices conjectured in (3.12) are to be completely consistent then the part 
arising from the fusion of S must provide some appropriate zeros in order to cancel the 
unwanted poles of the Toda theory factor. In the next section we show how this can occur 
by solving the bootstrap equation for the a m _i and c m cases. The obstruction for dealing 
with the b m and d m algebras in the same way, apart from the problem that the vector 
particle is not elementary, is that the relevant spectral decompositions of the .R-matrices 
are not known. 

S-matrices associated to the vector representations of the classical Lie algebras have 
recently been proposed in [13]. However, this reference uses an incorrect form for the 
spectral decomposition (3.16) and misses the crucial factor of t in the above equations 
which render the result invalid for c m . It also asserts that the masses of the theories are 
those of the corresponding classical Toda theory; whereas our results show that they are 
actually those of the related classical Toda theories (3.20). 

4. The solution of the bootstrap for a m _i and c m 

In this section we will explicitly solve the bootstrap for the a m _i and c m theories. 
This involves the resolution of two interlocking problems. Firstly we must find the spectral 
decompositions of the .R-matrix on the fundamental representations and then we must show 
how the scalar factor provides the correct analytic structure on the physical strip. 

Fortunately, for the representations that are necessary to construct the ^-matrices, 
there is well-established technology for finding the spectral decompositions of the associ- 
ated .R-matrices. We explain how the spectral decompositions are found in appendix A. 
However, knowing the spectral decompositions is not enough for our purposes, they only 
tell us the form of the .R-matrix up to a scalar function of the spectral parameter. In fact, 
we must solve the bootstrap equations for which the overall scalar factors are determined 
and indeed crucial for the correct form of the S'-matrix. We derive these scalar factors in 
appendix B. 

Our results are the following. We denote by R ah (u) the .R-matrix between the a th 
and 6 th fundamental representations (where the highest weight of the a th fundamental 
representation is u a = e± + e2 + • • • + e a ), and we choose without loss of generality b > a. 
For a m _i we have 

min(m— b,a) 

R a \u) = Zf{u) (-)" + Vf(^ b+fe+ ^- fe , (4- 1 ) 

fc=0 
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with 

k min(m— b,a) 

P a k b (u) = f[{2p + b-a} J] {-2p-b + a}, 
p=i p=fc+i 

a b— 1 a 

Zf(u) = Hl[{2j+2k-a-b} H {- 2 p-b + a}. 

3=1 k=l p=min(m — b,a) + l 



where we take uj m = ujq = and 



For Cm we have 



Zf (u) = (^^) [} \\{2j + 2k-a-b}{a + b-2(m + l)-2j -2k} 
\smuj j =lk=1 



x n { 2 P i> ■ (,\. 

p=min(m — b,a) + l 



The full S'-matrix is then equal to 

jab 
»(fc) 



with 

a b 



(4.2) 



= sin(W2 + An)_ 
sin a; 



min(m — b,a) a—j 
j=0 fc=0 

where 

= II^P + 6 " a > II{ 2 ( m + 1) + 2? - a - 6} 

p=l q = l 

min(m— b,a) a 

x ]"| {-2p-b + a} Yl {-2(m + l) -2q + a + b}. 

p=j + l q = k + l 

ab a 6—1 



(4.5) 



Sf'(u)=Y ab (u)R a \u\ (4.6) 



Y ab (u) = f[ f[ Y (u + l(2j + 2fc - a - b - 2)). (4.7) 
i=ifc=i V 19 J 

where Y(u) is given by (3.6) for a m _i and (3.9) for c m . In (4.6) we have explicitly indicated 
the dependence on the parameter k via uj in (2.15). The normalization factors Z\ h {u) and 
Z% b (u) arise from solving the bootstrap equations starting with R(u) and are crucial for 
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the correct form of the S-matrix since they can provide zeros on the physical strip and 
ensure that the unitarity condition is satisfied. To illustrate the latter point we notice that 

a b , 1 v 
R ab {u)R ba {-u) = \\ Y[q (u+— (2j + 2k-a-b-2) , (4.8) 
i=ifc=i V ^ / 

where g(u) is defined in (2.12). So using (4.7) along with (3.7) and (3.8) we deduce that 

S ab (u)S ba (-u) =I b ® I a . (4.9) 

We are now in a position to investigate the analytic structure of Sfy (u) on the physical 
strip. For a m _i there are no poles or zeros on the physical strip, whereas for c m there are 
no poles but if a + b > m there is a series of simple zeros at 

o + 6-2j+2 - = 12 +6 _ m _ (410) 

2(m+l) V ; 

A consistent S'-matrix is made by appending a suitable minimal Toda factor which 
provides the necessary pole structure. Generalizing the case for a m _i one is tempted to 
try the Gross-Neveu (GN) type Ansatz 

S$ ) (u) = X ab (u)S? b ) (u), (4.11) 

where X ab (u) is the minimal Toda factor of the associated algebra (3.19). However, we 
shall argue below that this is inconsistent for the non-simply-laced algebras due to the 
appearance of spurious poles on the physical strip which cannot be explained. Rather we 
shall find that a consistent S'-matrix is given by the tensor product form which generalizes 
the S-matrix of the Principal Chiral Model (PCM): 

S(t, W = X ab (u)S^(u) ® (4.12) 

For this S-matrix the particles transform in the reducible representations V LJa ® V LUa where 
V LUa is the representation of the algebra with highest weight u a . 

Let us verify the above statements for algebras a m _i and c m . The former case is easy 
to discuss. The minimal Toda factor is associated to the algebra a!^-i- 

a+b-l 

X a \u)= J] (j + 1)(j-1), (4.13) 

3 = |a-b| + l 
step 2 

with the notation 

7TU I 7T j \ 



U) = ) 9 ( , (4-14) 
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where tg = to in this case. The pole in X ab (u) at (a + b)/m (if a + b < to) or 2 — (a + b)/m 
(if a + 6 > to) corresponds to particle a + 6ora + 6 — to, respectively, in the direct channel. 
One can verify directly using (4.1) that the residues at these poles are P^ a+b and P Wo+b _ m , 
respectively, as required for consistency. The pole at \a — b\/m corresponds to the particle 
| a — b\ in the crossed channel. The element S ab (u) also exhibits double poles at 

a + b — 2k 

u= , k = 1,2, . . . ,min(a,b) - 1. (4.15) 

to 

These double poles are understood in terms of the Coleman-Thun mechanism [4] which 
was originally posited in order to explain the double poles in the sine-Gordon S-matrix. 
Essentially, the poles are Landau singularities of the ordinary sort, which appear as poles 
in two-dimensions and branch-points in four-dimensions, and are associated to 'on-shell' 
diagrams. In the present example, figure 1 gives the on-shell diagram causing the k th double 
pole in (4.15) for the case a, b < to/2 (other cases are found by crossing). We therefore 
conclude that both the generalized GN- and PCM-type S-matrices are consistent for a m _i. 




1. On-shell diagram giving double pole. 

(2) 

The scalar factor for c m is the d m+l minimal Toda S'-matrix. This S'-matrix can be 
understood as the S-matrix of the d^ +2 excluding the spinor and anti-spinor particles. So 

a+b-l 

X ab (u)= J] (j + l)(j-l)(2g-j + l)(2g-j-l), (4.16) 

j = |o-6| + l 
step 2 

where g = m + 1 is the dual Coxeter number of c m and the bracket notation is defined 
in (4.14). So X ab (u) in general exhibits poles of order 1,2,3 and 4 on the physical strip. 
Consider a GN-type S-matrix. The analytic structure of Sfy (u) is deduced by combining 
that of X ab (u) with the zeros (4.10) of the S ab (u) part. Rather than analyzing the S-matrix 
in full detail we shall exhibit a singularity which has no explanation in terms of the particle 
spectrum. To this end one notices that X a ' m+1 ~ a (u) has a double pole at u = 1/2 which is 
explained in terms of the on-shell diagrams in figure 2 involving the spinor and anti-spinor 
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particles of the d^ +2 theory. 7 S ab (u) has a simple zero at u = 1/2; hence Sfy(u) has a 
simple pole at u = 1/2 which cannot be explained in terms of the particle spectrum. At 
this stage we are forced to conclude that there cannot be a consistent S-matrix of the 
GN-type. 8 



m+l-a 




m+l-a 

2. On-shell diagram giving double pole in X a,m+1 ~ a (u) for <4n+2- 

The PCM-type of .S-matrix, however, does not have this spurious pole since there is 
an extra zero to cancel it. The analytic structure on the physical strip is illustrated in 
figure 3 where the crosses represent poles. 



a+b<m+l 



XXX ■ 




XXX xxx ■ 

1/2 


■ xxx 

1 


a+b-m+1 






XXX 


::»<!*» 







1/2 


1 


a+b>m+l 






xxx : 


: : xxx : f xxx : 





1/2 1 



3. Poles in the PCM-type S- matrix for c m . 

Notice that there are only simple poles and double poles. Some of the simple poles 
can be understood in terms of direct or crossed channel resonances of the particles and we 
consider these ones first. It is useful to distinguish three cases. 

7 I am grateful to Patrick Dorey for explaining this point. 

8 It might be possible to enlarge the set of particles and then close the bootstrap, however, 
the new particle would transform in a reducible representation involving non-fundamental 
representations. 
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(i) a + b < m + 1. In this case there are four simple poles. The simple poles at 
(a + b)/2(m + 1) and 1 — \a — b\/2(m + 1) correspond to the exchange of particles a + b and 
\a — b\ in the direct channel, respectively, with residues proportional to P^ a+b and P^ ^. 
The simple poles at \a — b\/2(m + 1) and 1 — (a + b)/2{m + 1) correspond to the exchange 
of | a — b\ and a + b in the crossed channel, respectively. 

(ii) a + b = m + 1. In this case there are two simple poles. The simple poles at 
\a — b\/2(m + 1) and 1 — \a — b\/2(m + 1) correspond to the exchange of particle \a — b\ in 
the crossed and direct channel, respectively. In the latter case the residue is proportional 
toP W|o _ b| . 

(iii) a + b > m+ 1. In this case there are four simple poles, however, only two of them 
correspond to direct or cross channel resonances. The simple poles at \a — b\/2(m + 1) and 
l — \a — b\/2(m + 1) correspond to the exchange of particle \a — b\ in the crossed and direct 
channel, respectively. In the latter case the residue is proportional to P u , _,, . 

We now turn to the remaining singularities. Firstly ifa + 6>m + l there are simple 
poles at (a + b)/2(m + 1) and 1 — (a + b)/2(m + 1). Secondly there are double poles are 
in the union of the set 

a + b — 2k 

u= ^Tf ^TV' fc = l,2,...,min(a,&)-l, (4.17) 

2[m + 1) 

and its crossed version, with the exceptions that fora + 6>m + lw = (a + b)/2(m + 1) 
and u = 1 — (a + b)/2(m + 1) are simple rather than double poles and when a + b = m + 1 
u = 1/2 is regular. 

There exists at least one on-shell diagram at the positions of the double poles, which 
involve particles in the spectrum and which lead to double poles on kinematical grounds. 
These on-shell diagrams are the same as those in figure 1 (plus the crossed versions). If 
a+b > m + 1 one would expect to see triple poles at (a+6)/2(m + l) and 1 — (a+b)/2(m+l) 
coming from the factor X ab (u), whereas in fact the S'-matrix only has simple poles. It is 
rather unconventional to have an S'-matrix, some of whose simple poles do not correspond 
to direct or cross channel resonances. However, [14] discusses a recent example of such an 
eventuality in non-simply-laced Toda theories. These strange occurrences are explained 
by a generalization of the Coleman-Thun mechanism [4], in the sense that there are on- 
shell diagrams which at first sight seem to lead to higher order poles but which on closer 
inspection actually only lead to simple poles because some sub-S'-matrix element in the 
diagram has a zero. We shall find that a similar mechanism is at work in the example 
under discussion, except that here the softening of the singularity is not due to some sub- 
S'-matrix element having a zero but rather is intimately bound up with the fact that the 
particles carry internal quantum numbers and one must sum all the quantum numbers of 
particles on internal lines. 
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The diagram which looks like it would lead to a triple pole at (a + b)/2(m + 1) (and 
1 — (a + b)/2(m + 1) by crossing) is illustrated in figure 4 (with the incoming particles 
coming from the left), where the grey blob is the S-matrix element s rn+1 ~ a ~ b,rn+1 ~ a ~ b (u) 
evaluated at u = uq = 1 — (a + b)/2(m + 1) and the black blobs represent projection 
operators. 




4. On-shell diagram for a + b > m + 1. 

As it stands, standard kinematical arguments would indicate that the diagram yields 
a double pole. However the S'-matrix element s m+1 ~ a ~ b,m+1 ~ a ~ b (u) has a simple pole 
at uq corresponding to the exchange of 2(m + 1) — a — b in the direct channel; hence the 
expected singularity is a triple pole. However, this argument is too naive and we must 
investigate the residue more closely. The behaviour of the internal S'-matrix element in 
the vicinity of the pole to 0(u — uq) is 



-<m+l — a — b,m+l— a — b 



(«) = 

(a^ + (u- uq) b ^^j + ( u ~ u o) b ^ F ^j ' (4 ' 18) 

where P' = ^uj 2{m+1) _ a _ b , a and b^ are constants and the sums are over the highest weights 
of the representations that appear in the tensor product V UJm+1 _ a _ b <E> V UJm+1 _ a _ b . The 
relevant observation is that the representation with highest weight uj2( m +i)-a-b does not 
actually appear in the tensor product V u!a <8> V UJb (remember that a + b > m + 1). So if we 
isolate the group-theoretic contributions from one of the factors in the tensor product to the 
S-matrix we conclude that the terms of the form PoP'P , where Po represents the product 
of the two projection operators on the external legs and hence is a homomorphism from 
Vu a ®Vu b to V Um+1 _ a _ b <S>V ulm+1 _ a _ b , must vanish. Hence when the diagram is evaluated one 
actually picks up the contribution from the central S-matrix element at 0(u — uo), rather 
than at O (u — Mo) -1 • I n other words the contributions from the more singular diagrams 
are zero after one has summed over all the intermediate states. So the actual singularity 
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is a simple pole rather than a triple pole since each triangle contributes a simple pole from 
its kinematical factors. 

It is important to be aware of the fact that the S-matrix X ab (u) also has fourth order 

(2) 

poles which can be explained solely in terms of the set of particles in the d m+1 spectrum. 
These singularities are reduced to poles of second order in the PCM-type S'-matrix. In 
order for this to occur the residues of the diagrams corresponding to the fourth order poles 
must vanish. Such an eventuality is possible because the particles carry internal quantum 
numbers — as we have seen from the preceding arguments. Unfortunately, an analysis 
of the residues is a rather formidable problem and is not within the scope of the present 
article. 

Before we move on to discuss the rational limits of the S-matrices we first pause to 
consider the d m and b m cases. As we have pointed out for these algebras the vector is not 
the elementary particle, rather it is the spinor (and anti-spinor) particle. Nevertheless our 
S-matrices should generate a subset of the full S'-matrix. One finds that in these cases the 
bound states are not associated to the fundamental representations, as in the a m _i and 
c m cases. In fact the a th particle transforms in the reducible representation [15] 

a-2j>0 

W a = V Ua _ v . (4.19) 
i=o 

Unfortunately the spectral decompositions have not been found for these representations. 
5. The rational limit 

The rational limit of the S-matrices are obtained by taking k — > oo and the resulting S- 
matrices are those of some well-known quantum field theories. Furthermore one can show 
by explicit computation that the 'important' analytic structure of the S-matrix is not 
affected by the limit, in the sense that no poles or zeros from the S factor can wander onto 
or off the physical strip. The rational S-matrices are actually invariant under the group 
associated to the Lie algebra in question (on the contrary the trigonometric S-matrices of 
the last section are invariant under the quantum group). 

The rational limit of the i?-matrices is easily obtained by taking the k — > oo limits of 
(4.1) and (4.4), which means replacing (4.3) with 

{x}=(x + tgu)/2, (5.1) 

and taking the P's to be the Lie algebra, rather than quantum group, invariant homomor- 
phisms. This latter limit is taken because as k — > oo the deformation parameter of the 
quantum group q — > — 1 and so the quantum group reduces to the Lie algebra. 
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The factor Y tgitk (u) also has a good limit: 



Y tg,oo{ u ) 



i r(i-f)r(^ + f) 
fcr(i + s)r(i + £-5 




(5.2) 



Using these results we find for vector-vector scattering 



S(oo)(w) 



Lli^lllii+l) 




p 2wi + 




(5.3) 



for a m _i and 




(5.4) 



for the other algebras. 



A careful comparison of the rational limits of our S-matrices (5.3) and (5.4) with those 
of the principal chiral model [3], 9 in which the particles transform in a tensor product of 
fundamental representations of the algebra, confirms that 



for a m _i and c m (and b m and d m for a = b = 1). This relation explains our use the 
nomenclature 'PCM- type' for the S-matrix of (4.12). 
For a m _i we can also consider the S-matrix 



which is the S-matrix of the SU(m) Gross-Neveu model [3,16]. 

The fact that the Gross-Neveu S'-matrix for the non-simply-laced algebras violates 
the bootstrap was noted in [3] . It has been claimed [3] that the S'-matrix of the principal 
chiral model for c m violates the bootstrap due to the appearance of simple poles with no 
explanation in terms of direct or cross channel poles. However, the analysis of the last 
section shows how this problem is resolved and the poles can be properly understood via 
a generalization of the Coleman-Thun mechanism. 

9 In comparing our expressions with this reference it is important to notice that the S'-matrices 
defined there are equal to ours up to a permutation of the outgoing particles. 



^cm(«) = S(t,oo) W = s ( t } H ® S&^X^u), 



(5.5) 



S a G b N (u) = S& ) (u) = S& ) (u)X ab (u), 



(5.6) 
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6. 



The situation at k = 1 



In this section we consider what happens to the S'-matrices when the level k = 1. 
It is known from the study of quantum groups that when q is a root of unity, so k is 
an integer, that the representation theory has to modified. This modification is imple- 
mented automatically by moving to the IRF picture and means that the usual rule for 
decomposing tensor products is truncated on the set of representations which correspond 
to highest weights of level k. This would entail appropriate modifications of the spectral 
decompositions. When k = 1 only level one representations survive. It is significant that 
in the cases where the particles are associated to higher level representations, i.e. for the 
b m and d m theories where there are representations at level 2, then the representations 
(4.19) are reducible and contain a level 1 component, so that all the particles remain in 
the spectrum. 10 

As we mentioned, the appropriate way to write down the S'-matrix for k = 1 is via 
the IRF formalism. At k = 1 the allowed weights are simply A*(l), that is 

{0,wi,w 2 ,...,w ra -i}, fora m _i, {0,^i,w m }, for b m , 

{0,wi,w 2 ,...,w ra }, for c m , {0, u>i, u>m-i, ^m}, for d m . 

For each of the representations W a (recall from (4.19) that they are in general reducible 
for b and d series) one can draw an admissibility diagram consisting of a set of nodes 
for each element of A*(l) with a joined to b by an oriented link if a — b is a weight of 
the representation. For example, in figures 5 and 6 we give the admissibility diagrams 
for the vector spinor and anti-spinor representations of d m and for the vector and spinor 
representations of c m 

One notices immediately a qualitative difference between the admissibility diagrams 
of the simply-laced and non-simply-laced algebras. In the former case the admissibility 
diagrams are trivial in the sense that the space of paths on the diagram from a given 
starting point and given length is just one-dimensional. This means that the resulting 
solution of the YBE is trivial: 

Sg?)(«) = l. (6.2) 

This means that the GN-type S'-matrix for the simply-laced algebras is related to the 
PCM-type S-matrix: 

SUu) = SH !,(«), (6.3) 

10 This observation seems also to be true for the known solutions of the YBE for the exceptional 
algebras as well. 
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vector (1) 



spinor(m-l) anti-spinor (m) 




5. Admissibility diagrams for d m with k = 1. 



vector (1) spinor (m) 




6. Admissibility diagrams for c m and k = 1. 

and furthermore 

= S ( ?)(«) = (6.4) 

the minimal Toda S-matrix. 

On the contrary, for the non-simply-laced algebras the admissibility diagrams are not 
trivial when k = 1. In other words at = 1 the S'-matrix is still non-diagonal and a 
relations like (6.3) and (6.4) do not hold. 

7. Discussion 

We have constructed factorizable S-matrices for trigonometric solutions of the YBE 
for the vector representations of all the classical Lie algebras. By appending a suitable 
CDD factor the S-matrices have singularities corresponding to new states. For the a m _i 
and c m algebras, where the relevant spectral decompositions are known, one can solve the 
bootstrap equations to find that the spectrum consists of particles transforming in the 
fundamental representations of the algebra. For b m and d m the situation is less clear and 
in any case the vector particle is not the fundamental particle; rather one would expect 
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this to be the spinor (and anti-spinor for d m ). 

The simplest Ansatz for the S-matrix, generalizing that of the Gross-Neveu model, 
fails for c m due to the existence of singularities on the physical strip which cannot be 
explained in terms of the spectrum of particles. A generalized principal chiral model 
Ansatz, for which the particles transform in a tensor product of representations, is shown 
to be consistent and all singularities on the physical strip can be accounted for. 

It has been argued that the a m _i trigonometric S-matrices describe certain integrable 
deformations of some coset conformal field theories [6,8]. Consider the coset conformal field 
theory Xk x xi/x^+i, f° r some Lie algebra x. There is a particular relevant operator in the 
theory of holomorphic dimension 

A = l--=— ? , (7.1) 

where g is as before the dual Coxeter number of x, which leads to a massive theory 
with higher spin integrals of motion. The natural conjecture is that the S'-matrix of this 
integrable theory is precisely the generalized PCM-type S'-matrix S®£ ^(u) [17]. Notice 
that the x\ x Xk/xk+i cases are described by the GN-type S-matrix for the simply-laced 
algebras only (see section 6). The conjecture could perhaps be placed on a better footing 
by employing the technology of the Thermodynamic Bethe Ansatz to investigate the ultra- 
violet limit of the theories. The relation of the trigonometric solutions to the YBE and 
Bethe Ansatz systems has been considered in [3,18,19]. 

Finally, the a m _i trigonometric S- matrices have been proposed to describe the scat- 
tering of solitons in complex a^_ 1 Toda field theory [5,20]. This has been established via 
semi-classical techniques. It is now known that all complex affine Toda theories admit soli- 
ton solutions [21], and is natural to ask whether the trigonometric S-matrices for the other 
algebras describe the scattering of these solitons. This may be true for the simply-laced 
algebras, for which the soliton S-matrix would be the generalized GN-type S-matrix. For 
the non-simply-laced algebras the situation is much less clear. Experience with the real 
Toda theories in these cases shows that the resolution of the problems may be surprisingly 
subtle and perhaps will require new types of solution to the YBE appropriate to a set of 
particles whose mass ratios depend on a coupling constant. 

I would like to thank Patrick Dorey for many conversations on S-matrices and also 
Gustav Delius for some useful discussions. Also I would like to thank Niall Mackay for 
pointing out an error in an earlier version of the manuscript which allowed for a significant 
improvement. 
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Appendix A: The spectral decompositions for a m _i and c m 

In this appendix we derive the spectral decompositions of the .R-matrices on the 
fundamental representations of the a m _i and c m algebras. We follow the approach of [22] 
although our method is also a direct translation into the quantum group of [23] which 
found the spectral decompositions of the rational R- matrices for c m . 

The first point to make is an .R-matrix cannot be associated with any two representa- 
tions of the quantum group. It is necessary that the representations are affinizable in the 
language of [22]. It can be shown that all the fundamental representations of a m _i and 
c m have this property. However, it is not true for the fundamental representations of the 
other algebras, where the affinizable representations are reducible in general [10,15]. 

Now we apply the technology of [22,23] to find the spectral decompositions of the 
.R-matrices on the fundamental representations. One first constructs the Tensor Product 
Graph (TPG). This is a graph is constructed by letting the irreducible components of the 
tensor product ® V v be the nodes joined by a link if V\ and V a have opposite parity and 
V a C adjoint<g>Vx- The parity of an irreducible component V\ is defined to be ±1 according 
to whether V\ appears symmetrically or anti-symmetrically in the tensor product (in the 
limit q — > 1). 

For V (x>a <g> V UJb (with b > a without loss of generality) of a m _i, using the notation 
j, k = V Uj +uj k and j = V LJj we have the TPG 

a, b <-> a — 1,6+1 •••<-> a — min(m — 6, a), b + min(m — 6, a) . (A.l) 

For c m we have for a + b < m 

a, b <-> a — 1,6+1 

I I 
a — 1,6— 1 <-> a — 2, 6 

i i 

l,6-a+ 1 <-> 6-a + 2 
I 

6 — a 

whilst if a + 6 > m then the graph truncates at the (m — 6 + l)th column. 
The spectral decomposition of the i?-matrix has the form 

R a \x) = Y,P^W^ (A.3) 



<-> l,a + b — 1 <-> a + 6 
I 

<-> a + 6-2 



(A.2) 
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where the sum is over the representations that appear in the tensor product V LOa ® V^ b and 
hence is a sum over nodes of the TPG. x is the multiplicative spectral parameter. If there 
is an arrow from v to p on the TPG then the coefficients p^(x) and p v {x) satisfy 

p M (x) _ xqW 2 -x- x q I{y)l2 
Pu (x) ~ x"V (At)/2 -xgW 

where q is the deformation parameter of the quantum group and I(p) = (p + 2p) • p (where 
p is the sum of the fundamental weights) is the eigenvalue of the quadratic Casimir on the 
representation with highest-weight p. 

The spectral decompositions follow from applying the rule (A. 4) recursively from 
u a + Ub (for the case of c m it is important to notice that the result is independent of the 
path). With 

x = exp(zA-u), q = — exp(— iu), (A. 5) 

one finds (4.1) and (4.4) up to an overall multiplicative factor. 

Appendix B: Zeros and the bootstrap 

In this appendix we explain how the solution of the bootstrap equations leads to the 
zeros in (4.1) and (4.4). From the spectral decompositions of the .R-matrices we know that 

5^« 6 )ocP Wc , (B.l) 

for some appropriate value of u c ab , in which case we say that there is a fusion ab c. This 
also implies that there are fusions ac — > b and be —> a where a is the charge conjugate 
particle (for a m _i a = m — a whilst for the other algebras a = a). The identity 

Kb + + <-c = 2, (B.2) 

holds between the fusing parameters. The bootstrap equations give the S'-matrix elements 
of c in terms of a and b which follows from the fusion relation between the .R-matrices: 

R dc (u) = (l b ® R da (u + u\S) {R db {u - < c -) <g) 7°) , (B.3) 

restricted on the left and right to the subspace V UJc C V Ub ® V Ua . In the above u = 1 — u 
and I a is the identity on V LOa . Consider the case d = a evaluated at u = — u b a5 then 

R ac (-< 5 ) = {i b ® ^ aa (o)) (R ab (-u c ab ) ® I a ) (B.4) 
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(A.4) 



using (B.2). But from the spectral decompositions (4.1) and (4.4) one finds _R aa (0) oc 
I a ® I a and that R ab (— u c ah ) has no component proportional to P Wc , hence 

R ac (-<- c ) = 0. (B.5) 

We now apply (B.5) and (B.3) recursively to find zeros of R ah {u). Starting from R 11 ^) 
((3.13) and (3.16)) (B.5) implies that R 12 (-l/tg) = 0. For all the algebras except a m _i 
crossing symmetry would imply R 12 (l + 1/tg) = in addition. By the recursive use of 
(B.3) we find that R ab (u) for b > a has a set of zeros at 

u= — (a + b-2j - 2k), j = l,2,...,a, k = 1, 2, . . . , b - 1, (B.6) 
tg 

and for all the algebras except a m _i their crossed values as well. 

The remaining zeros in (4.1) and (4.4) are accounted for in a different way. When 
a + b > m then we saw in appendix A that the TPG gets truncated which implies an 
additional set of zeros at 

u= —(2k + b-a), k = m-b+l,m-b + 2,...,a. (B.7) 
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